Initial eccentricity and constituent quark number scaling of elliptic flow in ideal and 

viscous dynamics 
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In the Israel-Stewart's theory of dissipative hydrodynamics, we study the scaling properties elliptic 
flow in Au+Au collisions. Initial energy density of the fluid was fixed to reproduce STAR data on 
cj> meson multiplicity in 0-5% Au+Au collisions, such that irrespective of fluid viscosity, entropy at 
the freeze-out is similar in ideal or in viscous evolution. Initial eccentricity or constituent quark 
number scaling is only approximate in ideal or minimally viscous (n/s — l/47r) fluid. Eccentricity 
scaling become nearly exact in more viscous fluid (n/s >0.12). However, in more viscous fluid, 
constituent quark number scaled elliptic flow for mesons and baryons split into separate scaling 
functions. Simulated flows also do not exhibit 'universal scaling' i.e. elliptic flow scaled by the 
constituent quark number and charged particles V2 is not a single function of transverse kinetic 
energy scaled by the quark number. From a study of violation of universal scaling, we obtain an 
estimate of QGP viscosity, r//s — 0.12 ± 0.03. The error is statistical only. Systematic error in rj/s 
could be as large. 

PACS numbers: 47.75.+f, 25.75.-q, 25.75.Ld 



I. INTRODUCTION 

Recently STAR collaboration published their analysis 
of elliptic flow (v 2 ) of identified particles in Au+Au col- 
lisions at v^— 200 GeV [![. In the relatively low px re- 
gion, pr < 2GeV, v 2 of identified particles, in each cen- 
trality bins studied, scales with transverse kinetic energy 
KEt — tut ~~ m - However, they did not find initial or 
participant eccentricity scaling of V2{pt) or v 2 (mx — m). 
At higher px, v 2 scales with constituent quark number. 
They also observed 'universal scaling', v 2 /(n q < v 2 > c fc)> 
elliptic flow scaled by the constituent quark numbers and 
charged particles v 2 , of different particle species in differ- 
ent collision centrality, scales with KEx/n q , the trans- 
verse kinetic energy per constituent quark number. Sim- 
ilar results are obtained in PHENIX measurements 
The constituent quark number scaling of elliptic flow sug- 
gests existence of an initial collective partonic state. Con- 
stituent quark number scaling of elliptic flow is naturally 
explained in a coalescence model 3, 4] , where elliptic flow 
of constituent quarks adds up. 

Elliptic flow in heavy ion collisions is best understood 
in a hydrodynamic model Elliptic flow measure 

the momentum anisotropy. In non-zero impact param- 
eter collisions, the reaction zone is spatially asymmetric. 
Differential pressure gradient convert the spatial asym- 
metry to momentum asymmetry. Approximate initial 
eccentricity scaling of elliptic flow is observed in ideal 
hydrodynamic model simulations for Au+Au collisions. 
Constituent quark number scaling, however is not indi- 
cated in ideal hydrodynamic simulation of Au+Au col- 
lisions. In the present paper, we investigate the scaling 
properties of elliptic flow in ideal and viscous hydrody- 
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namics. From the study of scaling violation, we also 
obtain a phenomenological estimate of QGP viscosity, 
i]/s = 0.12 ± 0.03. In recent years, there is considerable 
interest in viscosity of strongly interacting Quark-Gluon 
Plasma. Theoretical estimate cover a wide range, rj/s ~ 

0. 0-1.0 . String theory based models (ADS/CFT) give a 
lower bound on viscosity of any matter rj/s > 1/4tt 

In a perturbative QCD, Arnold et al Ml estimated rj/s ~ 

1. In a SU(3) gauge theory, Meyer [8| gave the upper 
bound rj/s <1.0, and his best estimate is ry/s=0. 134(33) 
at T = 1.165T C . At RHIC region, Nakamura and Sakai 
Q estimated the viscosity of a hot gluon gas as rj/s=0.1- 
0.4. Attempts have been made to estimate QGP viscos- 
ity directly from experimental data. Gavin and Abdel- 
Aziz [lj| proposed to measure viscosity from transverse 
momentum fluctuations. From the existing data on 
Au+Au collisions, they estimated that QGP viscosity as 
ry/s=0.08-0.30. Experimental data on elliptic flow has 
also been used to estimate QGP viscosity. Elliptic flow 
scales with eccentricity. Departure from the scaling can 
be understood as due to off-equilibrium effect and utilized 
to estimate viscosity [ll[ as, ?7/s=0.11-0.19. Experimen- 
tal observation that elliptic flow scales with transverse 
kinetic energy is also used to estimate QGP viscosity, 
rj/s ~ 0.09 ± 0.015 [3, a value close to the ADS/CFT 
bound. From heavy quark energy loss, PHENIX collab- 
oration [l3| estimated QGP viscosity rj/s ~ 0.1-0.16. In 
[iH [Hj], from analysis of RHIC data, an upper bound 
to viscosity is given rj/s < 0.5 [H, [l5j]. Recently in [l6| . 
from a detail analysis of centrality dependence of 4> meson 
multiplicity, mean px and integrated v 2 , QGP viscosity 
was estimated as, rj/s = 0.07 + 0.03 + 0.14, the first error 
is statistical, the second one is systematic. The large sys- 
tematic error is due to uncertain knowledge about various 
parameters of a hydrodynamics model, e.g. initial time, 
energy density profile, freeze-out condition,, finite accu- 
racy of computation etc. However, estimates of QGP 
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viscosity from experimental data must be treated with 
caution. Recently, in [llj , it was shown that viscosity to 
entropy ratio (rj/s) in Au+Au collisions depend on the 
collision centrality. While in central collisions, charged 
particle elliptic flow demand nearly perfect fluid, more 
viscous fluid is demanded in peripheral collisions. In the 
collision centrality, 0-60%, rj/s can vary from 0-0.17. 

The paper is organized as follows: in section|TTl we de- 
scribe the model: hydrodynamic equations for space-time 
evolution of the fluid, equation of state and initial condi- 
tions. Various scaling properties of elliptic flow in ideal 
and viscous fluid evolution are studied in section lHll In 
section llVl from violation of universal scaling, we obtain 
an estimate of viscosity to entropy ratio. Results are 
summarized in sectionfVl 



II. HYDRODYNAMIC EQUATIONS, 
EQUATION OF STATE AND INITIAL 
CONDITIONS 
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FIG. 1: Black circles are lattice simulation [23J for entropy 
density. The black line is the model EOS, obtained by para- 
metric representation to the lattice simulations and a hadronic 
resonance gas at low temperature. 



Israel-Stewarts theory of 2nd order dissipative rela- 
tivistic hydrodynamics is well established Briefly, 
in Israel-Stewart's theory, the thermodynamic space is 
extended to include the dissipative flows, relaxation 
equations for which are solved simultaneously with the 
energy-momentum, baryon number conservation equa- 
tions. More detail exposition can be found in (l9l. [2pj. 

Presently, we assume that in % /s ArAr =200 GeV, Au+Au 
collisions at RHIC, a baryon free fluid is formed. Only 
dissipative effect we consider is the shear viscosity, heat 
conduction and bulk viscosity is neglected. The space- 
time evolution of the fluid is obtained by solving, 

= 0, (1) 
Dtt^ = - — (tt'"' - 2? 7 V </ V /> ) 

TV 

- [u»tt" X +u v t: vX ]Du x . (2) 

EqUJ is the conservation equation for the energy- 
momentum tensor, T^ v = (e + p)u^u v — pg^ v + , 
£, p and u being the energy density, pressure and fluid 
velocity respectively, is the shear stress tensor. Eqf2] 
is the relaxation equation for the shear stress tensor tt^ v . 
In Eq[2J D — u^d^ is the convective time derivative, 
V«V> = ifV^H" + Vu^) - \(d.u)(g' iv - u^u") is a 
symmetric traceless tensor, ry is the shear viscosity and 
tv is the relaxation time. It may be mentioned that in 
a conformally symmetric fluid relaxation equation can 
contain additional terms (2lj . 

Assuming boost-invariance, EqsfT] and [5] are solved in 
(r = Vt 2 — z 2 ,x, y, r) s = 5 m fzf) coordinates, with the 
code " 'AZHYDRO-KOLKATA" \ developed at the Cy- 
clotron Centre, Kolkata. Details of the code can be found 
in Within 10% or less, AZHYDRO-KOLKATA sim- 
ulation reproduces Song and Heinz 's [2l| result for tem- 
poral evolution of momentum anisotropy e p . 



Eqs lll2l are closed with an equation of state p = 
p(e). Lattice simulations [22], [23[ indicate that the 
confinement-deconfinement transition is a cross over, 
rather than a 1st or 2nd order phase transition. In FigfTJ 
a recent lattice simulation [23| for the entropy densit y is 
shown. We complement the lattice simulated EOS |23j 
by a hadronic resonance gas (HRG) EOS comprising all 
the resonances below mass 2.5 GeV. In FigfTJ the solid 
line is the entropy density of the '"lattice +HRG'" EOS. 
The entropy density is obtained as, 

s = 0.5[1 - tanh(x)]s HRG + 0.5[1 - tanh(x)]si a ttice (3) 

with x = T ^ a ■ In the present simulation, we have used 
cross over temperature, T c =196 MeV and AT = 0.1T C . 
Compared to lattice simulation, entropy density in HRG 
drops slowly at low temperature. It is consistent with ob- 
servation in [23| , that at low temperature, trace anomaly, 
S y 4 p drops faster in lattice simulation than in a HRG 
model. It is difficult to resolve whether the discrepancy 
is due to failure of HRG model at lower temperature or 
due to the difficulty in resolving low energy hadron spec- 
trum on rather coarse lattice [23j. 

Solution of partial differential equations (Eqs lll2[) re- 
quires initial conditions, e.g. transverse profile of the 
energy density (e(x,y)), fluid velocity (v x (x,y),v y (x,y)) 
and shear stress tensor (7r M "(x, y)) at the initial time Tj. 
One also need to specify the viscosity (77) and the relax- 
ation time (r w ). A freeze-out prescription is also needed 
to convert the information about fluid energy density 
and velocity to particle spectra and compare with ex- 
periment. 

We initialize the fluid as in fig. We assume that the 
fluid is thermalized at Ti—Q.6 fm [5j and the initial fluid 
velocity is zero, v x {x,y) = v y {x,y) = 0. Initial energy 
density is assumed to be distributed as @ 
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FIG. 2: charged particle's integrated elliptic flow scaled by 
the participant eccentricity (i^/e) as a function of collision 
centrality. The lines with the filled circles, squares, triangles 
and inverted triangles are hydrodynamic predicted V2 in fluid 
evolution with rj/s=0, 0.08, 0.12 and 0.16 respectively. 



FIG. 3: charged particle's integrated elliptic flow scaled by 
the participant eccentricity (v2/e) as a function of g d ch , 
in ideal and viscous evolution. The blank symbols are the 
scaled flow when central energy density is fixed at eo=35.5 
GeV/fm 3 . The filled symbols are obtained in evolution when 
central energy density is reduced to keep multiplicity fixed. 
For ease of visibility, blank circles are shifted horizontally. 



e(b, x, y) = £j[0.75iVp art (b, x, y) + 0.25N coU (b, x, y)], 

(4) 

where b is the impact parameter of the collision. N part 
and N co u are the transverse profile of the average num- 
ber of participants and average number collisions re- 
spectively. Si is a parameter which does not depend 
on the impact parameter of the collision. The shear 
stress tensor was initialized with boost-invariant value, 
7r xx — n vv = 2r)/3Ti, tt xv =Q. For the relaxation time, 
we used the Boltzmann estimate t„ = 3r)/2p. Hydrody- 
namics also require a freeze-out condition. In the present 
work, we have assumed baryon free fluid. The chemical 
potential is zero throughout the evolution. We assume 
that the kinetic freeze-out occur at a fixed temperature 
TV =150 MeV (24|. Details of particle production in vis- 
cous hydrodynamics can be found in [2(|. Briefly, in- 
variant particle distribution from the freeze-out hyper- 
surface S M is obtained using Cooper-Frye prescription 

E 7tr = l^-= I d ^f^p) ( 5 ) 

dtp dyd^pT Jt, 

where f(x,p) is the one-body distribution function with 
a non-equilibrium correction, f(x,p) — f eq (x,p)[l + 

<t>{x,p)], 4>{x,p) = 2 (e+p)T* K^P^Pv 

We have simulated Au+Au collisions for four values of 
viscosity, (i) t]/s=0 (ideal fluid), (ii) r]/s = l/Ait « 0.08 
(ADS/CFT lower limit of viscosity) , (hi) r;/s=0.12 and 
(iv) rj/s—0. 16. In viscous fluid evolution, entropy is gen- 
erated. To obtain a common ground to compare between 
ideal and viscous evolution, we fix the initial energy den- 
sity of the fluid such that entropy at the freeze-out, both 



TABLE I: Initial central energy density (gj) and temperature 
(Ti) of the fluid in b=0 Au+Au collisions, for different values 
of viscosity to entropy ratio (rj/s). The predicted (j> meson 
multiplicity and mean pr are also noted. They should be 
compared with STAR measurements, = 7.95 ± 0.74 

ex 

and < v % > ex = 0.977 ± 0.064. 



T)/S 





0.08 


0.12 


0.16 


ei{GeV/fm 3 ) 


35.5 
± 5.0 


29.1 
± 3.6 


25.6 
± 4.0 


20.8 
± 2.7 


Ti (MeV) 


377.0 
±13.7 


359.1 
±11.5 


348.0 
±14.3 


330.5 
±11.3 


dNf 
dy 


7.96 


8.01 


8.22 


8.13 


< pT<j> > 


1.019 


1.062 


1.111 


1.174 



in ideal and viscous fluid evolution, is the same. This 
is done by reproducing the STAR data [26j], on meson 
multiplicity in 0-5% Au±Au collisions. In table 1, initial 
central energy density (ej) and temperature (Tj) required 
to fit STAR data on cf> meson multiplicity in 0-5% Au±Au 
collisions are noted. The error in or in Ti corresponds 
to statistical and systematic uncertainty in STAR mea- 
surements [26] . The predicted <f> meson multiplicities and 
mean px are also shown in tablc|U Irrespective of fluid 
viscosity, initial energy density (or temperature) of the 
fluid can be tuned to reproduce STAR measurements for 

J TV T <p 

(j) meson multiplicity, = 7.95 ± 0.74 (statistical 

and systematic error included). As entropy production 
increases with viscosity, with increasing viscosity, fluid 
require less and less initial energy density. For example, 
compared to ideal fluid in minimally viscous fluid, ini- 
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FIG. 4: (color online) Elliptic flow scaled by the participant 
eccentricity («2/e) as a function of transverse momentum. 
The different symbols are for tt~ , K + , p, (j> and f2 in evo- 
lution of QGP fluid with viscosity to entropy ratio (a)r?/s=0, 
(b)77/s=0.08, (c) ?7/s=0.12 and (d) r]/s=0.16. 



tial energy density is reduced by ~ 18%. In fluid with 
viscosity 77/5=0. 16, reduction is even more, ~40%. In ta- 
bleU we have also noted the hydrodynamic predictions 
for <p meson mean px ■ < Pt > increases in viscous evolu- 
tion and one find that predictions overestimate the STAR 
measurement, < p T >= 0.977 ± 0.064 (statistical and 
systematic error included) for rj/s >0.12. Experimen- 
tal data on (f> meson multiplicity and mean px in 0-5% 
centrality Au+Au collisions are simultaneously explained 
only for rj/s <0.12. 



III. SCALING PROPERTIES OF ELLIPTIC 
FLOW 

A. Initial eccentricity scaling 

We first investigate the initial eccentricity scaling of 
integrated elliptic flow in ideal and viscous dynamics. It 
is now realized that the experimentally measured scaling 
of integrated V2 with multiplicity or with collision cen- 
trality is not in agreement with hydrodynamics. While 
hydrodynamics predicts approximate scaling (2lj . in ex- 
periments scaling is violated [13, H^|. As discussed in 
[29l | , violation of the scaling can be understood as an in- 
dication of incomplete thermalization. However, there is 
a serious flaw in the arguments of [2§| . In [231 j a relation 
between inverse of the Knudsen parameter i and mul- 
tiplicity per unit the transverse area ^ ^ , was obtained 
under the assumption that the total particle number is 
conserved throughout the evolution. The assumption was 
justified with the observation that particle density is pro- 
portional to entropy density and entropy is conserved. In 
a viscous evolution, entropy is generated and initial and 
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FIG. 5: (color online) same as in Fig(3j but as a function of 
transverse kinetic energy. 



final state entropy are not same and the assumption is 
clearly violated. Only in systems with very small vis- 
cosity, the assumption may be valid, but not in systems 
where sufficient entropy is generated. Explicit numerical 
simulations show (see table 1) that entropy generation 
can be substantial, e.g. ~ 20%, 30% and 50% in fluid 
evolution with viscosity to entropy ratio, ?7/s=0.08, 0.12 
and 0.16. 

In FigJU present model predictions for the eccentric- 
ity scaled charged particle's elliptic flow, in ideal and 
viscous dynamics, as a function of participant numbers 
(Np ar t) are shown. Except for the very peripheral colli- 
sions, V2/S, in ideal or viscous dynamics, approximately 
scales with collision centrality. In very peripheral colli- 
sions, N part <50, V2 in viscous hydrodynamics is more 
than that in ideal hydrodynamics. However, applica- 
bility of hydrodynamics in very peripheral collisions is 
questionable [5j. Apparently, the results ap pea r to con- 
tradict simulation studies of Song and Heinz [2l| . In [2l[ , 
Song and Heinz studied initial eccentricity scaling of inte- 
grated V2 as a function of charged particle multiplicity per 
unit transverse area, They obtained approximate 

scaling both in ideal and viscous fluid, however, scal- 
ing violation is somewhat larger in viscous dynamics. In 
|2l| fluid was not initialized to have similar multiplicity. 
Viscosity reduces elliptic flow. However, since entropy 
is generated in viscous evolution, particle multiplicity is 
increased. A part of larger scaling violation in viscous 
dynamics is due to enhanced multiplicity and reduced 
elliptic flow. In the present simulation, fluid was initial- 
ized to have similar multiplicity, both in ideal and viscous 
dynamics. Since energy density is reduced, viscous sup- 
pression is also comparatively less. Consequently, scaling 
violation is not large in viscous fluid evolution. This is 
explicitly shown in Fig[31 In Fig[31 the blank circles, 
squares, up triangles and down triangles corresponds to 
eccentricity scaled elliptic flow in hydrodynamic simu- 
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FIG. 6: (color online) Elliptic flow scaled by the constituent 
quark number for ir~ , K + , p, <j> and in 20-30% Au+Au 
collision as a function of transverse momentum scaled by the 
constituent quark number (pr/n q ). Scaled flow in deal and 
viscous t?/s=0.08, 0.12 and 0.16 fluid are shown separately. 



lation with viscosity to entropy ratio ri/s=0, 0.08, 0.12 
and 0.16 respectively. The initial central energy density 
is fixed at £o=35.5 GeV/ fm 3 . With increasing viscos- 
ity, increases but elliptic flow decreases. The filled 
symbols in Fig |3] are the scaled flow when multiplicity is 
constrained by reducing central energy density in more 
viscous fluid (see tabled . With the exception of very pe- 
ripheral collisions, dispersion between the filled symbols 
is comparatively less than that between the blank sym- 
bols. For a given viscosity, elliptic flow is less suppressed 
when multiplicity is constrained. 

To investigate the initial eccentricity scaling of differ- 
ential elliptic flow of different particle species in different 
collision centrality, we have computed elliptic flow for 
7i- , K+, proton, 4> and fl in 0-10%, 10-20%, 20-30%, 30- 
40%, 40-50% and 50-60% Au+Au collisions. Note the 
particle list contain non-strange and strange mesons and 
baryons. In a hydrodynamic model, elliptic flow of iden- 
tified particles, in different collision centrality, approx- 
imately scales. To show the scaling between different 
particle species, in FigPJJ scaled elliptic flow for (i) 7r~ in 
0-10%, (ii) K+ in 10-20%, (hi) in 20-30%, (iv) p in 30- 
40% and (vi) f2 in 50-60% Au+Au collisions are shown. 
Only a few collision centralities are chosen to show the 
scaling behavior clearly. Results are essentially same if 
some other centralities are chosen. In the present work, 
we have neglected resonance decays. Resonance contri- 
bution reduces elliptic flow, mostly at low pt 30]. Pions 
are most affected by resonance decays. In the low px 
range, < pr < 1 GeV, V2 is reduced by ~ 0-30%. At 
Pr > 1 GeV, resonance contribution to V2 is negligible 
[30j | . Mass ordering of the flow is clearly evident at low 
Pt, more V2 for lighter particles. But at large pr mass 
effect is reduced. One observe that as a function of pt 
elliptic flow do not show initial eccentricity scaling. Same 
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FIG. 7: (color online) same as in Fig|6]but as a function of 
constituent number scaled transverse kinetic energy KET/n q . 



results are shown in Fig {5] as a fu nction of transverse ki- 
netic energy, KEt — \J m 2 + p T — m. As a function of 
KEt, initial eccentricity scaling of elliptic flow is only 
approximate in ideal fluid evolution. Scaling gets better 
with viscosity. For example, in ideal hydrodynamics, at 
KE T ~ 2 GeV, v 2 /e for fi in 50-60% Au+Au collisions 
is ^45% larger than the same for 7r~ in 0-10% Au+Au 
collisions. The difference is less in viscous evolution, e.g. 
fl flow exceed that of tt by, - 30%, ~15% and -7% 
for ?7/s=0.08, 0.12 and 0.16 respectively Simulations in- 
dicate that initial eccentricity scaling of elliptic flow is 
better observed in viscous fluid evolution than in ideal 
fluid evolution. The result is important as it indirectly 
indicate that in Au+Au collisions, a low viscous fluid is 
produced. Otherwise, in experiment more exact initial 
eccentricity scaling would have been observed. 



B. Constituent quark number scaling 



While constituent quark number scaling is natural in 
a coalescence model, [H, IH, it is not indicated in a hy- 
drodynamic model. In a hydrodynamic approach, quark 
degrees of freedom are not manifestly evident. However, 
since in experiment, constituent quark number scaling is 
observed, and if hydrodynamics do faithfully represent 
the experimental data, one would expect, to be consis- 
tent with experiment, meson and baryon flow in the ratio 
2:3 . 

In FiglHl elliptic flow for 7r~, K + , p, <fi and ft, scaled 
by the constituent quarks number in 20-30% Au+Au col- 
lisions are shown as a function of pr /n q . Quark number 
scaling of elliptic flow, as a function of pT/n q naturally 
occur within a coalescence model Q. However, experi- 
mentally scaling is better observed as a function of trans- 
verse kinetic energy. In FigjTl constituent quark number 
scaling as a function of KET/n q are shown. Hydrody- 
namics predictions for constituent quark number scaling, 
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FIG. 8: (color online) initial eccentricity and constituent 
quark number scaled elliptic flow (v2/n q e) for ir~ , K + , p, <f> 
and Q as a function of KEt /n q , in QGP fluid evolution with 
(a) r)/s=0, (b) r]/s=0.08, (c) tj/s=0.12 and (d) rj/s=0.16. 
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FIG. 9: (color online) charged particle «2 and constituent 
quark number scaled elliptic flow (v2/n q e) for ir~ , K + , p, (/> 
and f2 as a function of KEr/n q , in QGP fluid evolution with 
(a) r?/s=0, (b) ?7/s=0.08, (c) r?/s=0.12 and (d) r?/s=0.16. 
Approximate 'experimental' scaling function, as obtained in 
PHENIX experiment is shown as the solid line. 



as a function of pt or KEt are similar, but the scaling 
seems to work better as a function of KEr/riq. At a fixed 
centrality, either as a function of pt or KEt, constituent 
quark number scaling is only approximate in ideal fluid 
dynamics. In a viscous fluid, constituent quark number 
scaling seems to work better but separately for mesons 
and baryons. For fluid viscosity ?7/s=0.12-0.16, the n q - 
scaled v 2 of baryons and meson is approximately in the 
ratio 2:3 at KET/n q > 1 GeV. The result is expected. 
As shown in FigGH for ?y/s=0.12-0.16, elliptic flow for 
mesons and baryons are approximately same at large pt ■ 
When scaled by the constituent quark number, naturally, 
the baryon and meson flow are in the ratio 2:3. 

In FigEl elliptic flow scaled by the initial eccentricity 
and constituent quark numbers as a function of trans- 
verse kinetic energy are shown. As before 7r~, K + , 
p, cj) and fTs are from 0-10%, 10-20%, 20-30%, 30-40% 
and 50-60% centrality Au+Au collisions. ^ scaling is 
also an approximate in ideal fluid, or minimally viscous 
(?//s=0.08) fluid evolution. For example, at KEx/uq ~ 
.6 GeV, V2/n q e of ir~ in 0-10% collision and that of ft in 
50-60% collision differ by ~30% in ideal and by -10% 
in minimally viscous fluid evolution. At larger KEt 
the difference is even more. For higher values of viscos- 
ity ry/s=0.12 and 0.16, at KE T /n q > 0.5 GeV, v 2 /n q e 
for mesons and baryons scale separately. The result is 
consistent with nearly perfect eccentricity scaling in vis- 
cous (t7/s=0.12-0.16) fluid evolution (see Fig(5|) and sep- 
arate constituent quark number scaling for meson and 
baryons (see Fig(7j) . Violation of initial eccentricity and 
constituent quark number scaling in ideal or viscous hy- 
drodynamics is also consistent with STAR and PHENIX 
experiments. Experimental data do not show the scaling. 



C. Universal scaling of elliptic flow 

Both the STAR QJ and PHENIX collaboration ob- 
served universal scaling of elliptic flow. Elliptic flow 
scaled by the constituent quark number and the charged 
particle v 2 , for different particle species in different col- 
lision centrality is approximately a single function of 
the quark number scaled transverse kinetic energy. The 
STAR collaboration observed the scaling for K® A, S in 
0-10%, 10-40%, 40-80% Au+Au collisions (see Fig. 12d 
in Q). In PHENIX collaboration showed that both 
in Au+Au and Cu+Cu collisions, the universal scaling is 
approximately valid for a large number of particle species 
in a variety of collision centrality (see Fig. 3 in 0). 

In Fig. IH1 we have examined the universal scaling in 
the hydrodynamic model. Hydrodynamical simulations 
for V2/{n q * 3.2* < v 2 > c h), for 7T, K, p, tfi and fl, as a 
function of KET/n q are shown. < v 2 > c h is the appro- 
priately weighted pt integrated flow for ir, K and p. In 
experiments, charged particles list is mainly populated 
by these three species. In Fig|9l the scaled flow, for tt, 
K, p, <fi and f2 are shown. Contrary to experiments, in 
hydrodynamic simulations, irrespective of fluid viscosity, 
elliptic flow does not follow the universal scaling. 



IV. UNIVERSAL SCALING VIOLATION AND 
QGP VISCOSITY 

In previous sections we have studied several scaling 
properties of elliptic flow in ideal and viscous dynamics. 
It appears that ideal or viscous hydrodynamic simula- 
tions for elliptic flow in Au+Au collisions do not show 
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the scaling properties exhibited in experiment. However, 
until now we have not compared hydrodynamic simula- 
tions for elliptic flow with experiments. In this section, 
we try to obtain an estimate of QGP viscosity by con- 
fronting hydrodynamic simulations for elliptic flow with 
experiment. 



TABLE II: Departure of hydrodynamic simulations for elliptic 
flow of 7r, K , p, 4> and Q, in the collision centrality 0-60% from 
the PHENIX universal scaling function. 



rj/s 


TT 


K 


P 





n 





0.87 


2.20 


3.04 


5.36 


6.17 


0.08 


0.25 


0.67 


0.76 


1.53 


1.62 


0.12 


0.03 


0.09 


0.36 


0.51 


0.67 


0.16 


0.05 


0.09 


0.91 


0.97 


1.50 



As stated earlier, PHENIX collaboration observed uni- 
versal scaling of elliptic flow for different particle species, 
in different centrality ranges of Au+Au and Cu+Cu col- 
lisions In Fig|9l the solid line represents the approx- 
imate scaling function as obtained in the PHENIX ex- 
periment 0. Experimental scaling function does have 
fluctuations, more at large KET/n q , which is presently 
ignored. Comparing the experimental scaling function 
with hydrodynamic simulations, we observe that, (i) in 
ideal fluid simulations, elliptic flow is over predicted for 
all the particle species, (ii) elliptic flow is also over pre- 
dicted in minimally viscous fluid evolution, though the 
difference between theory and experiment is less than 
that in ideal fluid evolution, (iii) flow is under-estimated 
in evolution with 7y/s=0.16, and (iv) simulated flows show 
least departure from the experimental flow for 77/5=0. 12. 

To measure the departure of simulated flow from the 
experimental 'universal scaling function', we define a 
scaling violation function F. 

F = ZM h s C ale d )*-V2scale d } 2 (6) 

where v^scaied IS approximate scaling function as ob- 
tained in PHENIX experiment (the solid line shown in 
Figinj) an d (vacated)* is the hydrodynamical simulations 
for the scaled flow v^lijiq * 3.2* < v 2 > c h), in the i — th 
centrality bin. As stated earlier, we have neglected reso- 
nance contribution. For low mass particles, present simu- 
lations for V2 in low pt, Pt _1 GeV, may not be reliable. 
Accordingly, we compute F in the pr region 1 < px < 3 
GeV. The quantity F is computed separately for parti- 
cle species tt, K, p, <j> and fi, in 0-10%, 10-20%, 20-30%, 
30-40%, 40-50% and 50-60% Au+Au collision. The re- 
sults are shown in tablc|H] Departure from universal 
scaling function is mass dependent. Hydrodynamic sim- 
ulated flow for heavy mass particles deviate more from 
the universal scaling function than that of a lighter parti- 
cle. Scaling violation is also viscosity dependent. Scaling 
violation is more in ideal fluid evolution than in viscous 
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1 \ 

0.00 0.05 0.10 0.15 0.20 
x\/s 

FIG. 10: x 2 / N f or the combined data set, tt + + tt~, K + +K~ 
and p + p as a function of viscosity to entropy ratio r)/s. The 
solid line is a parabolic fit to the x 1 /N values. 

fluid evolution. It may also be noted that for the viscosity 
to entropy ratio, rj/s— 0.12, the violation is minimum. 

In Fig llOl the scaling violation function F, summed 
over all the particles (tt, K , p, (f> and Q) is shown as a 
function of rj/s. The solid line in Fig|10] is a parabolic 
fit to F. From the minimum of the fit, we estimate vis- 
cosity to entropy ratio, for which simulated flows show 
least departure from the PHENIX experimental scaling 
function as rj/s = 0.12 ± 0.03. 

The present estimate rj/s = 0.12 ± 0.03 is well within 
the upper bound of viscosity obtained in [3, HH • Within 
the error, the estimate also agrees with other estimates 
obtained in a hydrodynamical model [ill [l2|. However, 
one must treat the phenomenological estimates with cau- 
tion. As noted in [16j , systematic error in hydrodynamic 
evaluation of viscosity could be large. Indeed, in 
from a simultaneous fit to <fi meson multiplicity, mean 
Pt and integrated v 2 , viscosity to entropy ratio was esti- 
mated as, rj/s = 0.07 ±0.03 ±0.14, where the 1st error is 
statistical and the 2nd one is the systematic error. Large 
(~ 200%) systematic error arises due to uncertainty in 
initial time, initial energy density profile, initial fluid ve- 
locity, freeze-out condition, finite accuracy of computer 
code etc. Even then as noted in [16| the source of sys- 
tematic error is not exhaustive. 

Before we summarise our results, few comments are 
in order. Present simulations indicate that elliptic flow 
scaling is better observed in viscous hydrodynamics than 
in ideal hydrodynamics. The result is contradictory to 
our expectation. In viscous hydrodynamics, Cooper-Frye 
formula for invariant distribution ( , d Jf — ), in addition 
to the thermal part, includes a viscous correction [20|. 
Since viscous correction introduces a microscopic scale, 
it is expected that any scaling property that would have 
been observed in ideal fluid, will worsens. Simulation 
studies indicate that the contribution of the viscous part 
to elliptic flow can be large ~ 50% [H, H3|. However, 
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viscous contribution to elliptic flow does depend on vari- 
ous parameters, e.g. if freeze-out temperature is lowered 
viscous correction decreases [20]. Viscous flux (n t21 ') de- 
creases rapidly with time and its contribution to Cooper- 
Frye decreases if fluid freezes out at later time or at 
lower temperature. Unlike in [2(], [lH, |33| , in the present 
simulations, a lattice based EOS with a confinement- 
deconfinement cross-over transition is used. There is no 
mixed phase and viscous fluxes decreases more rapidly 
than in 1st order phase transition. We have checked that 
in our simulations, at the late stage, fluid behaves nearly 
as an ideal fluid and viscous correction in Cooper-Frye 
is very small. Viscosity however changes the freeze-out 
surface. Changed freeze-out surface, but small viscous 
correction, may possibly be the reason for obtaining bet- 
ter scaling in viscous dynamics than in ideal dynamics. 

V. SUMMARY 

To summarize, in ideal and viscous hydrodynamics, we 
have studied scaling properties of elliptic flow in Au+Au 



collisions. Fluid was initialized to reproduce experimen- 
tal <f> meson multiplicity in 0-5% centrality Au+Au col- 
lisions. Elliptic flow in ideal or minimally viscous fluid 
evolution, show only approximate initial eccentricity scal- 
ing as a function of transverse kinetic energy. Initial ec- 
centricity scaling is better observed in more viscous fluid 
evolution. Constituent quark number scaling is also an 
approximate scaling in ideal and minimally viscous dy- 
namics. In more viscous fluid, constituent quark number 
scaling seems of be obeyed separately for mesons and 
baryons. Elliptic flow in ideal or viscous dynamic also 
do not show experimentally observed 'universal scaling'. 
From a study of universal scaling violation in ideal and 
viscous evolution, we also obtain an estimate of QGP vis- 
cosity, r]/s = 0.12 ± 0.03. However, the error does not 
include systematic error due to uncertain initial condi- 
tions in a hydrodynamic evolution. Systematic error in 
7]/s could be as large as ~ 200%. 
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